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Abstract 



"O - 

^ . First-principles calculations within the weighted density approximation (WDA) were performed 

o, 

CJ , for ground state properties of ferroelectric perovskites PbTiOa, BaTiOs, SrTiOa, KNbOa and 

^ ' KTaOs. We used the plane- wave pseudopotential method, a pair distribution function G based on 

> ■ 

^S] ' the uniform electron gas, and shell partitioning. Comparing with the local density approximation 

a^. 

CD . (LDA) and the general gradient approximation (GGA), we found that the WDA significantly im- 

^D I proves the equilibrium volume of these materials in cubic symmetry over both the LDA and GGA; 

^^ ' Ferroelectric instabilities calculated by the WDA agree with the LDA and GGA very well; At the 



experimental ferroelectric lattice, optimized atom positions by the WDA are in good agreement 
with measured data; However the WDA overestimates the strain of tetragonal PbTiOs at exper- 
imental volume; The WDA overestimates the volume of fully relaxed structures, but the GGA 
results are even worse. Some calculations were also done with other models for G. It is found that 



K> ! a G with longer range behavior yields improved relaxed structures. Possible avenues for improving 

C^ ' the WDA are discussed. 



I. INTRODUCTION 

Since the early 1980s, first-principles calculations based on the density functional theory 
(DFT) have been implemented to compute diverse properties of piezoelectric and ferroelec- 
tric materials. The main difficulty within the DFT is how to treat the exchange-correlation 
(xc) energy accurately, because the exact form of it remains unknown. The local density ap- 
proximation (LDA), in which the xc energy density (exc) depends only on local charge density, 
dominated these calculations due to its simplicity and surprising success. When performed 
at the experimental volume, the LDA predicts many properties of ferroelectric materials, 
such as phonon frequencies, ferroelectric phase trausi.ious. polarizatiou. elasfcity. etc.. with 
extraordinary accuracy JL,]- However it is well known that the LDA overestimates the bind- 
ing energy, and underestimates the bond length by 1-2%, which results in the calculated 
equilibrium volume normally 3-6% less than experiment. On the other hand, ferroelectric 
properties are extremely sensitive to volume. For example, the ferroelectric instabilities in 
BaTiOa J2| and KNbOa pl, m \^ m] are severely reduced, if not totally eliminated, at the 
LDA zero pressure volume. Even at the experimental volume, the LDA incorrectly predict 
certain properties, e.g., it overestimates the strain of the tetragonal PbTiOi ll| and the 
static dielectric coustaut .. QQ. aud it uudcrcsti.ates the baud gap flQ. 



liQ, 



which includes the density gra- 



The generalized gradient approximation (GGA) 
dients, is the natural next step beyond the LDA. Generally speaking, the semi-local GGA 
tends to improve upon the LDA in many aspects, especially for atomic energies and struc- 
tural energy differences Ij, |l3|, [l^. For ferroelectric properties and band gaps, the GGA 
normally predicts very similar results to the LDA at experimental volume. However the GGA 
tends to overestimate the bond length by about 1% p, |l5[. To illustrate this we present a 
full relaxation of tetragonal PAmm PbTiOs within both the LDA (Hedin-Lundqvist [16[) 



and the GGA (The PBE 12| version was used throughout this paper.) using the linearized 
augmented planewave method with local orbital extensions (LAPW+LO) [l2|. As seen in 
TableHl at the experimental volume, the GGA predicts a strain better than the LDA, but the 
fully relaxed GGA structure is much worse than that of the LDA. We also used the ABINIT 



package [18|, |19|, which is based on pseudopotentials and a planewave basis. For the fully 
relaxed tetragonal PbTiOs, the LDA predicts a volume of 60.34 A'^, and c/a = 1.042, while 
the GGA gives a volume of 70.54 A^, and c/a = 1.24. The ABINIT results are in excellent 



agreement with the LAPW results. The failures of the LDA and GGA indicate that more 
complicated non-local approximations are needed. 



The weighted density approximation (WDA) J20|, |2l|, |22], within which e^c depends on 
charge density over a finite region, was advanced in the late 1970's. The WDA assumes 
that any inhomogeneous electron gas can be regarded as continuously being deformed from 
a homogeneous electron gas, so the real pair-distribution function of the inhomogeneous gas 
can be replaced by that of a homogeneous gas at every point with certain weighted density. 
By constructing a model xc hole, the sum rule is fulfilled. One problem that both the LDA 
and the GGA have is the self- interaction error, i.e., the imperfect cancellation of Hartree 
and xc terms in the one-electron limit. The self-interaction-correction (SIC) method 



predicts band gaps of transition-metal oxides in good agreement with experiments 2J] . It 
demonstrates the importance of SIC. However the SIC method is not efficient because its 
potential is orbital-dependent. The WDA is free of self-interaction for the one-electron limit, 
and its potential is directly a functional of the charge density, and not orbital dependent. 
For the non-local (uniform) limit, it reduces to the LDA. It has been reported that the WDA 



significantly improves atomic energies [22, |2^ and the equilibrium volume of some simple 



bulk solids |2^ over the LDA. It is promising to apply the WDA for ground state properties 
of more complicated ferroelectric perovskites. 

The selected perovskites PbTiOa, BaTiOs, SrTiOs, KNbOa and KTaOa are both techni- 
cally and theoretically important ferroelectric materials. They are also prototypes of more 
complicated and interesting ferroelectric relaxors such as (1— a;)Pb(Mgi/3Nb2/3)03-xPbTi03 
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. They have previously been 



27 



28|. In this paper, first 



(PMN-PT) and (1 - a;)Pb(Zni/3Nb2/3)03-xPbTi03 (PZN-P': 

studied extensively to understand ferroelectricity |2, S 14, S 

a brief overview of the WDA formalism will be presented in section II. Then we will report 

the WDA results of ground state properties of these perovskites in section III. In section 

IV, we discuss avenues for improving the WDA focussing on the symmetry under particle 

exchange in the pair correlation function. 



II. FORMALISM 

The general form of the xc energy in the DFT scheme can be expressed as (All equations 
are in atomic units, h = m = e"^ = 1.) 

E^n] = I J n{r)dr J ^0^dr\ (1) 

here the xc hole density nxc(i", r') is defined as 

n.,(r,r') = n(r') A^^,(r, r'; A) - l]dA 
Jo 

^n(r')K,(r,r')-l], (2) 

where f7xc(^' ^') i^ ^^^ coupling-constant averaged pair-distribution function of a system with 
a density n(r) 29|. This function is symmetric 

9:,{r,r')=g:^{r',r), (3) 

and the xc hole satisfies the following sum rule 

J n,,{r,r')dr' = -1. (4) 

Although the function g^c of a uniform gas is known with high accuracy based on Monte Carlo 
simulations [2l|, the exact form of g^^{r,r') is still elusive for an arbitrary inhomogeneous 
system. Different density approximations make different approximations to g^c- In the WDA 
approach (^"^.(r, r') is approximated by a model function G (e.g., a uniform type), 

G[|r-r'U(r)]=^,",(r,r')-l, (5) 

where the parameter n(r) is the weighted density, and it can be determined from the sum 
rule 

fn{r')G[\r-r'\,n{r)]dr' = -l. (6) 

By definition, the function G is not symmetric 

G(r,r')^G(r',r). (7) 

Since n(r) is only a function of r, not r', it leads to wrong asymptotic behavior at the low 



density limit 



which we will discuss in section IV. 



The natural choice of G would be obtained from g for a uniform electron gas, but since the 
systems of interest are inhomogeneous, other ty pes of model function G could be better J3l|, 



22] • Among them, the Gunnarsson- Jones (GJ [3l|) ansatz ensures the correct behavior at 



GGJ(r,n) = Ci(n){l-exp(-[— -]-5)}, (8) 

C2[n) 



large distance 

while Rushton et al. (RTZ J33|) used a simple Gaussian function 

G^^^{r,n)=c^{n)eM-[^f). (9) 

which resembles the exact one for the short distance limit, and it is also a good approximation 
to the uniform G. The parameters Ci and C2 can be determined from the following conditions 

n f G{r,n)d^r = -1, (10) 



n 



G (r, rij ^3 _ uqjjj 



r 



d'r = e^rin), (H) 



where e^°^{n) is the xc energy density of a homogeneous gas with density n. 
The WDA xc energy is 

ET>] = ljJ^^^^G[\r-r'\Mr)]drdr', (12) 

and the corresponding xc potential fxc(r) is the functional derivative of E^c, 

^r^(r) = -f;^ = Mr) + V2{t) + vs{r), (13) 



1 f n(r ) ^^1^ ^,| _^_^i j_, _ WDA/ 



where 

^i(r) = -y ^^G[|r-r'|,n(r)]rfr' = er^(r), (14) 

Mr) = ll^^^G[\r-r'\,n{r')]dr', (15) 

Examination of the above suggests that implementation of E^^^ and f^°^(r) could be 
cumbersome. However, in a plane-wave representation, by using the convolution theorem, 
these terms can be evaluated efficiently, as detailed in Ref. J26[. 

One subtle issue in the WDA implementation is shell partitioning J22]. The WDA scheme 
does not describe the exchange interaction between core and valence states very well because 



of the use of local n(r) in the model function G. Physically the range of integration of G is 
similar to the size of atoms, and there is no distinction between core and valence electrons to 
the contribution of E^J^^. As a result core and valence electrons dynamically screen valence 
electrons equally, which leads to non-zero exchange energy between core and valence states 
outside core regions. On the other hand, the LDA can give correct inter-shell contributions, 
since the LDA depends only on local density, and the core density vanishes outside core 



mm, 



regions. Based on this observation, a shell partitioning approach was proposed 

in which the valence-valence interactions are treated with the WDA, while core-core and 

core-valence interactions with the LDA. The total xc energy is written 

^xcN = E'^^^] + EZ^'^A - i^x'^^h.], (17) 



/' 



where n^ is the valence density and n is the total density. The sum rule becomes 

\n,{v')G[\v - r'|,n(r)] + n,{r)G[\r - r'|, n(r)]}cir' = -1, (18) 

where Uc is the core density. Ref 3J| shows why this sum rule must be used instead of 
a simpler one /{nv(r')G[|r — r'|, n(r)] = —1. The corresponding fxc for core and valence 
states will be different since E^c depends explicitly on both core and valence density. The 
detailed derivation and formulas can be found in the appendix of Ref. |3,5l]. We can still 
use the Hellmann-Feynmann theorem [;36] to determine atomic forces, and we found that 
the directly calculated forces agree with the numerical results of finite energy difference very 
well with the WDA shell partitioning. 

III. RESULTS 



A. Technical details 



y withi 



The WDA was implemented J26[ within a plane-wave basis pseudopotential method. The 
pseudopotentials are of the hard TrouUier- Mart ins type 37| . The semi-core states of metal 
ions include 3s and 3p states of K and Ti, 4s and 4j9 states of Nb and Sr, 5s and 5p states of 
Ta and Ba, and bd states of Pb. 2s states of O are also treated as semi-core. We used the 
WDA for valence states, the LDA for semi-core states, and pseudized lower states. Plane- 
wave basis sets with cut-off of 132 Ry were tested and found to be highly converged. A 
4x4x4 k-point mesh was exploited except for frozen-phonon soft-mode calculations of 



rhombohedral BaTiOs and KNbOa, where a denser 6x6x6 k-mesh was used because the 
ferroelectric double well depth is only about 2-3 niRy. The standard 3-point interpolation 
was employed to obtain n(r) with a logarithmic grid of increment nj+i = 1.25nj. We tested 
the convergence by the 6-point interpolation which makes negligible difference. For the 



function G, we used a uniform type 



38| for the following calculations. In the last part of 



this section we will also present results with other types of G. 

B. Equilibrium volume of cubic structures 

We first calculated the lattice constant of these materials constrained in cubic symmetry. 
The experimental values in Table |n] are those extrapolated to T = 0. The zero point 
corrections are expected to be less than 0.3%. As mentioned before, the LDA lattice constant 
is about 1-2% less than experiment, and this small error is big enough to make many 
ferroelectric properties incorrect. On the other hand, the GGA results are better than the 
LDA, but they are a little too large for BaTiOa, SrTiOa, and KTaOs. The WDA dramatically 
improves the lattice constant over the LDA, and it is also better than the GGA. Actually 
all these WDA lattice constants are very close to experimental data except for PbTiOa. At 
low temperature PbTiOs is tetragonal with a big strain of 6%, and other four perovskites 
are very close to the cubic structure. The extrapolation of the high temperature cubic data 
of PbTiOs therefore is not expected to be as reasonable as the other cases. 

C. Ferroelectric instabilities 

We studied the ferroelectric instability in tetragonal FbTiOs and rhombohedral BaTiOa 
and KNbOa. We displaced atom positions according to the experimental soft-mode distor- 
tion patterns at the experimental lattice. For tetragonal FbTiOa c/a = 1.063, while for 
rhombohedral BaTiOs and KNbOs, we neglected the tiny lattice distortion from the cubic 
structure. As mentioned before, the LDA and GGA describes the ferroelectric instability 
very well at the experimental structure, so one may hope that the WDA retains this good 
feature. Fig. ^ shows the calculated curves of energy versus ferroelectric displacement along 
[001] for PbTiOs and [111] for BaTiOg and KNbOs. In all cases, the WDA curves match 
with the LDA ones very well. It shows that the WDA can predict ferroelctric phase transi- 



tions as well as the LDA. On the other hand, the GGA predicts a smaller energy difference 
between paraelectric and ferroelectric states of PbTiOs than the LDA, but agrees well with 
the LDA for BaTiOg and KNbOs. 

D. Structural optimization 

Since there is no difficulty calculating atomic forces with the WDA approach, we per- 
formed structural optimization for tetragonal P4mm PbTiOa and rhombohedral P3ml 
BaTiOs and KNbOa. First we only optimized internal parameters at experimental lattices, 
and the relaxed atomic positions are shown in Tables ITTTl IIVI and IVl for PbTiOs, BaTiOa 
and KNbOs respectively. In all cases, the LDA, GGA, and WDA predict very similar results 
which are in good agreement with experiment except for KNbOs whose theoretical displace- 
ments from ideal positions are about 30% less than measured data [22|- The consistency of 
theoretical results suggests re-examination of these experimental data. 

We also did full relaxations for the above materials. At the experimental volume, the 
WDA predicts c/a = 1.106 for PbTiOs, which is close to the LDA value of 1.112, while 
the GGA gives c/a = 1.068. Compared with the measured value of 1.063, the WDA is 
not as good as the GGA in this case. For the fully relaxed structure, the WDA predicts 
a large c/a = 1.19 and a volume 8.2% bigger than experiment, which are poorer than the 
LDA, but still a little better than the GGA. To demonstrate that this is not because of 
poor pseudopotentials or other problems in the planewave method, we compared the energy 
difference between the LDA equilibrium structure and the WDA equilibrium structure. For 
the LDA, the planewave code predicts -11.8 mRy, while the LAPW gives -12.3 mRy. It proves 
that the planewave method is reliable. On the other hand, the WDA energy difference is 9.4 
mRy, which is similar in magnitude with opposite sign. We also calculated the equilibrium 
structures of rhombohedral BaTiOs and KNbOs, as shown in Table lVTl We can conclude that 
the LDA underestimates the fully relaxed volume by almost the same amount as constrained 
in the cubic structure, while the WDA overestimates it, and the GGA overestimates it even 
more. 



E. Results with other G functions 

Because there is no reason for the uniform G to be the best choice, we also tried other 
types of function G. We denote the uniform type as (a), the GJ type as (b), the simple 
Gaussian (RTZ) type as (c), and type (d) as G{r,n) = ci(n)exp(— [^7^]^). To clearly 
visualize them, we show these functions together at n = 1.0 in Fig. |21 The simple Gaussian 
type is a good approximation to the uniform G, the type (d) has longer interaction range 
than others, and the GJ type agrees with the uniform one when r > 1.745. 

We recalculated the equilibrium lattice constants of the cubic structure with the other 
three Gs. As shown in Table IVIH the newly calculated lattice constants are very similar. 
Compared with the previous (uniform G) results, the GJ type predicts slightly larger lattice 
constants, while the simple Gaussian form and type (d) predict slightly smaller ones. 

However the choice of G is more sensitive for the fully relaxed structure. Comparing with 
the uniform G, We found that the GJ type predicts even worse results, the simple Gaussian 
gives very similar results, and the type (d) is particularly better. For tetragonal PbTiOs at 
experimental volume, its optimized c/a becomes 1.078, and for the fully relaxed structure 
c/a is 1.092 and the volume is only 0.8% larger than experiment. However, as listed in Table 
IVIIIl the volumes of rhombohedral BaTiOs and KNbOa do not improve as much as that of 
tetragonal PbTiOa, they are still 2.4% and 1.3% larger than experiment respectively. 

We could further improve the WDA results by tuning the shape of function G, but it is 
hard to justify, and also it is difficult to make a single G to fit all properties of all materials. 
We need a simple and reasonable G, such as the uniform one, which can predict good ground 
state properties for both the cubic and the fully relaxed structures. Actually since different 
Gs predict very similar equilibrium volumes for the cubic structure, we expect their results 
for the fully relaxed structure should also be similar to the cubic structure, just like the 
LDA. So we focus on the WDA method itself. If one draws energy versus volume curves 
for the cubic and optimized tetragonal PbTiOs as shown in Fig. IHl one can see that the 
WDA curves for the cubic structure have similar curvatures to the LDA one, while for the 
tetragonal structure the WDA curves are too flat on the right side in comparison with the 
LDA. It means that for the relaxed structure the WDA predicts an energy increase smaller 
than expected on the large volume side. This suggests an asymptotic problem, as is discussed 
in the next section. 



IV. DISCUSSIONS AND PROSPECTS 

A. Asymptotic behavior 

An accurate functional approximation should fulfill at least some of the following con- 
ditions: (1) sum rule of the xc hole; (2) for a slowly varying density, it should recover the 
uniform gas limit; (3) absence of self-interaction; (4) asymptotic properties of xc energy and 
potential. One may verify that from the exact DFT expression (Eq. ^ the xc energy density 
exc(i") far away from the nucleus (low density limit) are 

lim e,e(r) ^ - — . (19) 

r—i-oo 2r 

Since g^^i^, r') is symmetric, one can derive that 

v^cir) = 2e^c{r)+V3{r), (20) 

and because ^"wrf vanishes exponentially in the above limit, 

lim v^,{r) ^ — . (21) 

The asymptotic conditions El and |^ are the consequence of requirement of cancellation 
of self-interaction in the Hartree terms. They look simple, but are very difficult to fulfill 
simultaneously. One can easily construct an energy or a potential satisfying the above 
conditions separately, but difficulty arises in making fxc(r) as the functional derivative of 
Exc[n{r)]. In the LDA and GGA, both exc(r) and fxc(r) go to zero exponentially for large r. 
Asymptotically they are less attractive than they should be. In the WDA, e^°^(r) satisfies 
condition Uni but the unsymmetrical function G (Eq. [7j) leads to fi(r) 7^ ^^2(1"), and f2(r) 
goes to zero exponentially also, giving that f!^°^(r) at large r to be the same as e^°^(r) 

lim t;WDA(r) ^ _1, (22) 

which is off by a factor of 1/2. 

The correct asymptotic behavior of e^°^(r) is due to the correct handling self-interaction 
in the WDA, while the wrong asymptotic behavior of w^^^(r) results from the violated 

symmetry of G under exchange (r, r') < > (r', r). In order to overcome this problem so that 

the WDA will behave closer to the exact DFT and eventually it will circumvent its failure 
for the fully relaxed structure, we can symmetrize G. 
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B. A new approach: symmetrization of the WDA 

A symmetric G satisfying Eq. El will make both e^°^(r) and v^J^^ijc) have the correct 
asymptotic behavior. The simplest approach is to make G depends on the weighted density 
at both r and r', for example, a sum form would be 

GSW°^(r,r') = ^{G[|r-r'U(r)] + G[|r-r'U(r')]}, (23) 

and 

<^SWDA(^^^,)^^SWDA(^/^^)^ (24) 

Now we have fi(r) = f2(r) = exc(r), and the conditions El and |^ are fulfilled naturally. 
The corresponding sum rule will be 

j n{r')G^'^^^{ry)dr' =]- j n{r'){G[\r - v'\,n{r)] + G[\r - r'ln{r')]]dv' = -I. (25) 

Self-consistent iteration can be used to determine n(r). Once nij) is known, E'^^^ is 
known also. The difficulty arises when one wants to calculate v^{y). If it can be determined 
efficiently, this symmetrized WDA will satisfy all the conditions listed in the beginning of 
last subsection, so this new WDA is very promising. It would overcome the problems that 
lead to shell partitioning because core and valence states would be distinguishable to screen 
valence electrons. 

V. CONCLUSIONS 

In conclusion, we calculated ground state properties of some common ferroelectric per- 
ovskites with the WDA. Compared with results of the LDA and GGA, the WDA describes 
properties of systems with the cubic symmetry or at experimental ferroelectric lattices very 
well, but it fails to predict good fully optimized structure. The symmetry problem of func- 
tion G could cause this failure. A new approach is proposed, and efforts must be taken to 
circumvent the mathematical difficulty. 
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TABLE I: Fully relaxed structures of tetragonal PbTiOa with the LDA and the GGA. Volumes 
are in A^, the numbers in parentheses are the deviations of strain and volume from experiment, 
and the experimental cja = 1.063. 



volume (Expt.) 



c/a 



volume (Relaxed) 



c/a 



LDA 
GGA 



63.28 
63.28 



L112 (+80%) 
L068 (+7%) 



60.36 (-4.6%) 
70.58 (+11%) 



1.051 (-20%) 
1.23 (+260%) 



TABLE II: Calculated LDA, GGA and WDA lattice constants in A for these ferroelectric materials 
in cubic state, compared with experimental data. Uniform electron gas G was used in the WDA. 
Numbers in parentheses are deviations from experiment. 



material 


LDA 


GGA 


WDA 


Expt. 


KNbOg 


3.960 (-1.4%) 


4.018 (+0.1%) 


4.011 (-0.1%) 


4.016 


KTaOg 


3.931 (-1.3%) 


4.032 (+1.2%) 


3.972 (-0.3%) 


3.983 


SrTiOg 


3.858 (-1.2%) 


3.935 (+0.8%) 


3.917 (+0.3%) 


3.905 


BaTiOs 


3.951 (-1.4%) 


4.023 (+0.6%) 


4.009 (+0.2%) 


4.000 


PbTiOg 


3.894 (-1.9%) 


3.971 (+0.1%) 


3.941 (-0.7%) 


3.969 



TABLE III: Optimized internal coordinates of tetragonal PbTiOa at the experimental volume 
(V = 63.28 A^) and strain {c/a = 1.063). Uz are given in terms of the lattice constant c. Uniform 
electron gas G was used in the WDA. 





LDA 


GGA 


WDA 


Expt. 


Uz(Ph) 


0.000 


0.000 


0.000 


0.000 


n.(Ti) 


0.539 


0.532 


0.539 


0.538 


n.(0i,02) 


0.615 


0.611 


0.614 


0.612 


n.(03) 


0.111 


0.105 


0.110 


0.117 
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TABLE IV: Optimized internal coordinates of rhombohedral BaTiOs for the experimental volume 
(V = 64.00 A^). Uz are given in terms of the lattice constant. Uniform electron gas G was used in 
the WDA. 





LDA 


GGA 


WDA 


Expt. 


n^(Ba) 


0.000 


0.000 


0.000 


0.000 


n.(Ti) 


0.488 


0.488 


0.489 


0.489 


n.(0i,02) 


0.511 


0.510 


0.509 


0.511 


n.(03) 


0.020 


0.018 


0.017 


0.018 



TABLE V: Optimized internal coordinates of rhombohedral KNbOa for the experimental volume 
{V = 64.77 A^). Uz are given in terms of the lattice constant. Uniform electron gas G was used in 
the WDA. 





LDA 


GGA 


WDA 


Expt. 


Uz{K) 


0.507 


0.509 


0.508 


0.5131 


n^(Nb) 


0.000 


0.000 


0.000 


0.0000 


n.(0i,02) 


0.019 


0.024 


0.017 


0.0313 


Uz{03) 


0.527 


0.522 


0.520 


0.5313 



TABLE VL Volumes (A^) and bulk modulus (GPa) of fully optimized tetragonal PbTiOa, rhom- 
bohedral BaTiOs and KNbOa. Uniform electron gas G was used in WDA. Numbers in parentheses 
are deviations from experiment. Experimental data are from Refs. l40l |. 







LDA 


GGA 


WDA 


Expt. 


PbTiOs 


V 


60.36 (-4.6%) 


70.58 (+11%) 


68.46 (+8.2%) 


63.28 




B 


83 


40 


42 


204 


BaTiOs 


V 


61.59 (-3.8%) 


67.47 (+5.4%) 


66.82 (+4.4%) 


64.00 




B 


148 


94 


99 


196, 156, 139 


KNbOa 


V 


61.96 (-4.3%) 


66.63 (+2.9%) 


66.04 (+2.0%) 


64.77 




B 


163 


98 


114 


138 
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TABLE VII: Calculated WDA (with four types of G denoted in subsection E) lattice constants 
in A for these ferroelectric materials in cubic state. Numbers in parentheses are deviations from 
experiment. 



material 


WDA (a) 


WDA (b) 


WDA (c) 


WDA (d) 


KNbOs 


4.011 (-0.1%) 


4.018 (+0.0%) 


4.006 (-0.2%) 


4.001 (-0.4%) 


KTaOg 


3.972 (-0.3%) 


3.979 (-0.1%) 


3.969 (-0.3%) 


3.964 (-0.5%) 


SrTiOa 


3.917 (+0.3%) 


3.924 (+0.5%) 


3.911 (+0.2%) 


3.908 (+0.1%) 


BaTiOg 


4.009 (+0.2%) 


4.015 (+0.4%) 


4.011 (+0.3%) 


4.004 (+0.1%) 


PbTiOg 


3.941 (-0.7%) 


3.948 (-0.5%) 


3.938 (-0.8%) 


3.932 (-0.9%) 



TABLE VIII: Volumes of fully optimized tetragonal PbTiOs, rhombohedral BaTiOs and KNbOs in 
A^ using the WDA with two types of G. Numbers in parentheses are deviations from experiment. 



WDA (a) 



WDA (d) 



PbTiOs 
BaTiOs 
KNbOa 



68.46 (+8.2%) 
66.82 (+4.4%) 
66.04 (+2.0%) 



63.79 (+0.8%) 
65.55 (+2.4%) 
65.60 (+1.3%) 
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FIG. 1: Total energy as a function of soft-mode displacement in tetragonal PbTiOs {c/a = 1.063), 
rhombohedral BaTiOs and KNbOs with the LDA, GGA, and WDA (uniform electron gas G). Here 
6 is the displacement relative to experiment. 
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FIG. 2: Four types of function G{r,n) with Vg = 1.0 (n = 0.2387). (a), (b), (c) and (d) correspond 
to the types defined in the last subsection of section III. 



18 



30 



25- 



,20- 






>,15 



m 



10 



5- 



1 r 


\ g 1 1 \ \ \ 1 1 TT 


1 r 


— 1 1 \ \ \ 


30 


- 


\ •/ 










> 






25 

















— Expt. 




20 




^-^/ 


«- 


— OLDA 






"^ 


\ /• 9 








w 


w »T t^r^ ytxf 


15 


\ 
\ 
\ 




o- 


-oWDA(d) 














~ \ 








10 


<> 


\ / yS^ 








^^ 


. ®^ ^*>x! "^w 








- 


V V«"''^^<f*v^. 






5 






*=- 


.^-— *^^1»^ 






1 1 1 1 ill 1 1 1 1 


'^r—*-'^ 


§ 


5 


60 65 „ 3 




70 7 




Volume (A ) 











FIG. 3: Energy as a function of volume for the cubic and the relaxed tetragonal PbTiOs, with 
the LDA, GGA, and WDA (two types of G). Dashed lines are cubic structures, and solid lines are 
tetragonal structures. 
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